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ABSTRACT 
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|"t | ■ We calculate the He I case B recombination cascade spectrum using improved 

q . radiative and collisional data. We present new emissivities over a range of elec- 

tron temperatures and densities. The differences between our results and the 
y—i . current standard are large enough to have a significant effect not only on the 

interpretation of observed spectra of a wide variety of objects but also on deter- 

Cs| . minations of the primordial helium abundance. 
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8 ■ 1. Introduction 
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Helium is the second most abundant element in the universe, and its emission and 
opacity help determine the structure of any interstellar cloud. Its abundance relative to 
hydrogen can be measured within a few percent since the emissivities of H I and He I lines 
have similar dependences on temperature and density. This makes it an indicator of both 
stellar and primordial nucleosynthesis (Pagel 1997). 



A good discussion of the history of calculations of the helium recombination spectra is 
given by Benjamin, Skillman, & Smits (1999, hereafter BSS99), who present new calculations 
- the current standard in the field. Yet much progress has been made since the work by Smits 
(1991, 1996) upon which the BSS99 results depend. We implement these improvements, 
present a new set of predictions, and compare our results with those of BSS99. The differences 
are large enough to impact continuing attempts to estimate the primordial helium abundance 
(Peimbert et al 2002). 
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2. The New Model Helium Atom 

The basic physical processes have been described by Brocklehurst (1972) and BSS99. 
Here we will describe the differences between BSS99 and our new numerical representation of 
the helium atom, which is a part of the spectral simulation code Cloudy (Ferland et al 1998). 
This model resolves all terms, nlS, up to an adjustable maximum principal quantum number 
n ma x, followed by a pseudolevel, n max + 1, in which all IS terms are assumed to be populated 
according to statistical weight and "collapsed" into one. We set recombinations into the 
collapsed level equal to the convergent sum of recombinations from n = n max + 1 to oo. In 
the low-density limit, the collapsed level increases the emissivities of our benchmark lines 
(the same 32 lines given in BSS99) by 0.4%, on average, with n max =100. The decays from 
states with I — n — 1 are most sensitive to this correction for system truncation. The strong 
optical line A5876 is corrected upward by 1.3%. At finite densities collisional processes force 
the populations of very highly excited states into local thermodynamic equilibrium (LTE). 
In this case the adequacy of the method used to compensate for truncation is unimportant. 
We find the corrections negligible for n e = 100 cm -3 and n max =100. Consequently, the 
uncertainties in the results presented in Section 3 are due to the uncertainties in atomic 
data, especially the often substantial uncertainties in collisional rates affecting terms not in 
LTE at given conditions. 

There are several differences in atomic data for radiative processes between BSS99 and 
the present work. The transition probabilities and radiative recombination coefficients are 
obtained from oscillator strengths and photoionization cross-sections. BSS99 uses the os- 
cillator strengths calculated by Kono & Hattori (1984). While these agree very well with 
the essentially exact oscillator strengths of Drake (1996), Drake presents a much larger set, 
up to and including n=10 and 1=7, which we adopt. Hummer & Storey (1998, hereafter 
HS98) have presented ab initio calculations of threshold photoionization cross-sections up to 
n=4. BSS99 uses cross-sections from TOPbase 1 (Cunto, 1993), while we use the more accu- 
rate HS98 values. The dominant remaining uncertainties in radiative data are in oscillator 
strengths involving low / states (with n > 10) and photoionization cross-sections for low / 
states (with n > 4). HS98 also illustrates the method, originally discussed by Seaton (1958), 
of calculating threshold photoionization cross-sections by extrapolating absorption oscillator 
strengths to the threshold energy of a given level. This method has been used in the present 
work, based on the oscillator strengths from Drake, to extend the ab initio cross-sections of 
HS98 to greater n. 

Differences in collisional data between BSS99 and the current work are also significant. 



1 http: / / vizicr.u-strasbg.fr/topbase/topbase.html 
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For low-n transitions for which there are ab initio calculations, BSS99 uses the collision 
strengths of Sawey & Berrington (1993). We replace these, where available, with the results 
of the close-coupling calculation by Bray et al (2000), which include continuum states not 
considered in the f?-matrix calculations by Sawey & Berrington. For /-changing collisions 
BSS99 uses two different treatments: Seaton (1962, hereafter S62) for low-/ transitions, and 
Pengelly & Seaton (1964, hereafter PS64) otherwise. Neither of these treatments allows 
for angular momentum transfers greater than one unit, and both apply when the projectile 
velocity is greater than the velocity of the bound electron. The r.m.s. electron and projectile 
velocities in conditions considered by BSS99, assuming proton colliders, are 



Zac , l3kT 

v e = and V proj = a , (1) 

n | m P 

where a is the fine structure constant, c is the speed of light, and Z is the screened nuclear 
charge. Dividing the latter by the former, we arrive at the expected value of the reduced 
velocity as a function of temperature and principal quantum number 



^ = Vprol = 71Q x 10 -5 nv /^|^ (2) 

For typical nebular temperatures this reduced velocity will be of order or less than unity 
for proton colliders for all n < 150; the treatment of PS64 is applicable only for greater n, 
and it greatly overestimates /-mixing cross-sections when used outside its range of validity 
(MacAdam, Rolfes, & Crosby 1981). 

Vrinceanu & Flannery (2001, hereafter VF01) give a (classical) theory of /-changing 
collisions and claim exact solutions in the limit that the intrashell transition is induced 
by slow distant collisions. Their treatment allows naturally for angular momentum changes 
greater than unity. (At a sufficiently high reduced velocity large angular momentum transfers 
are strongly suppressed and the theory goes to the optically allowed limit with which PS64 
is concerned.) We use equation 41 of Kazansky & Ostrovsky (1996) for the angle, A$, swept 
out by the projectile. A physical basis for the necessary large impact parameter cutoff in 
the theory follows from equating the Stark and quantum-defect precession frequencies. The 
Stark frequency is given by 



?>Zin 



(3) 



and the quantum-defect precession frequency (Hezel et al 1992) is given by 
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Fig. 1. — Percent difference between the emissivities calculated using two different Stark collision treat- 
ments, for several strong lines, as a function of n e . Left panel: a wide range of densities - as expected, 
there is no effect in either the low density extreme, because the collision rates are negligible, or the high 
density extreme, where the Stark collisions force the terms to LTE. The majority of lines arc most sensitive 
at densities found in stellar envelopes and quasar emission line regions. Right panel: the range of densities 
found in nebulae - several lines have a sensitivity to the Stark collision treatment of about 1%. 
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where 5i is the quantum defect, Z 1 is the charge of the projectile and b is the impact 
parameter. By setting u q d equal to u s , we obtain a maximum impact parameter, b max . The 
electron orbit precession will be faster than the Stark beating at larger impact parameters, 
so that transitions are increasingly less likely. To insure symmetry, we use the average u> q d 
of the initial and final levels. We use VF01 for /-changing collisions involving initial and 
final levels with I > 3, and like BSS99 we use the impact parameter treatment of S62 for 
/-change from s, p, and d levels. We use electron, proton and He + colliders for all transitions, 
taking n He + = 0.1 n p and n e = n p + n He +. Since S62, which describes electron collisions, 
is based upon the method of virtual quanta (see Jackson 1999), we can readily adapt it 
for the positive-ion collisions: The power spectrum of the time-dependent fields generated 
at the target atom by a passing charged projectile depends only on the projectile's charge 
magnitude, speed (not kinetic energy or mass separately) and impact parameter. The same 
considerations apply to PS64 and VF01 and have been implemented to allow for all three 
collider species. In calculating the necessary thermal averages we have assumed that the 
same temperature characterizes electrons, protons and He + ions. 

Figure 1 compares emissivities we predicted using the VF01 and PS64 theories. The 
predicted emissivities typically change by about 1% for nebular densities by using the theory 
of VF01 rather than that of PS64. The difference is much greater at high densities found, 
for example, in parts of quasars. 



In Figure 2 we compare our results with those of BSS99 for the case T e = 10 4 K and 
n e = 10 4 cm -3 . The average difference for the 32 emission lines is 4.6%. The greatest 
difference is for A4121, for which our emissivity is 25% greater. In general, agreement 
worsens with increasing density; at n e = 10 2 cm -3 , the average and greatest differences are 
1.6% and 6.4%, respectively, while at n e = 10 6 cm" 3 , we find differences of 7.0% and 35%. 
Agreement also worsens with increasing temperature. Table 1 presents emissivities for all 
of the temperatures and densities considered by BSS99. We believe that these results are a 
significant improvement. The application of these results to specific astrophysical problems 
will be the subject of future papers. 



3. 



Results 



We thank G. W. F. Drake for making available extensive tables of his calculations, D. 
Vrinceanu and M. J. Cavagnero for helpful discussions, and P. J. Storey, whose constructive 
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A comparison of the present results with those of BSS99 at T e — 10 4 K and n e = 10 4 cm~ 3 . 
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Table 1. He I Case B Emissivities. 
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Note. — Emissivities 4nj \/n e n He ^- are given in units 10 25 erg cm 3 sec 1 . The hydrogen density is 0.9 n e , and 
the helium abundance is one-tenth of the hydrogen abundance. 



